Analysis of dynamical tunnelling experiments with a Bose-Einstein condensate 
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Dynamical tunnelling is a quantum phenomenon where a classically forbidden process occurs, 
that is prohibited not by energy but by another constant of motion. The phenomenon of dynam- 
ical tunnelling has been recently observed in a sodium Bose-Einstein condensate. We present a 
detailed analysis of these experiments using numerical solutions of the three dimensional Gross- 
Pitaevskii equation and the corresponding Floquet theory. We explore the parameter dependency 
of the tunnelling oscillations and we move the quantum system towards the classical limit in the 
experimentally accessible regime. 
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I. INTRODUCTION 

Cold atoms provide a system which is particularly 
suited to study quantum nonlinear dynamics, quantum 
chaos and the quantum-classical borderland. On relevant 
timescales the effects of decoherence and dissipation are 
negligible. This allows us to study a Hamiltonian quan- 
tum system. Only recently dynamical tunnelling was ob- 
served in experiments with ultra-cold atoms PIQ- "Con- 
ventional" quantum tunnelling allows a particle to pass 
through a classical energy barrier. In contrast, in dynam- 
ical tunnelling a constant of motion other than energy 
classically forbids to access a different motional state. In 
our experiments atoms tunnelled back and forth between 
their initial oscillatory motion and the motion 180° out 
of phase. A related experiment was carried out by Steck, 
Oskay and Raizen in which atoms tunnelled from 

one uni-directional librational motion into another oppo- 
sitely directed motion. 

Luter and Reichl p| analyzed both experiments calcu- 
lating mean momentum expectations values and Floquet 
states for some of the parameter sets for which experi- 
ments were carried out and found good agreement with 
the observed tunnelling frequencies. Averbukh, Osovski 
and Moiseyev Q pointed out that it is possible to effec- 
tively control the tunnelling period by varying the effec- 
tive Planck's constant by only 10%. They showed one can 
observe both suppression due to the degeneracy of two 
Floquet states and enhancement due to the interaction 
with a third state in such a small interval. 

Here we present a detailed theoretical and numerical 



analysis of our experiments. We use numerical solutions 
of the Gross-Pitaevskii equation and Floquet theory to 
analyze the experiments and to investigate the relevant 
tunnelling dynamics. In particular we show how dynami- 
cal tunnelling can be understood in a two and three state 
framework using Floquet theory. We show that there is 
good agreement between experiments and both Gross- 
Pitaevskii evolution and Floquet theory. We examine 
the parameter sensitivity of the tunnelling period to un- 
derstand the underlying tunnelling mechanisms. We also 
discuss such concepts as chaos-assisted and resonance- 
assisted tunnelling in relation to our experimental re- 
sults. Finally predictions are made concerning what can 
happen when the quantum system is moved towards the 
classical limit. 

In our experiments a sodium Bose-Einstein conden- 
sate was adiabatically loaded into a far detuned optical 
standing wave. For a sufficient large detuning, sponta- 
neous emission can be neglected on the time scales of the 
experiments (160 fj,s). This also allows to consider the 
external degrees of freedom only. The dynamics perpen- 
dicular to the standing wave are not significant, therefore 
we are led to an effectively one-dimensional system. The 
one-dimensional system can be described in the corre- 
sponding two-dimensional phase space which is spanned 
by momentum and position coordinates along the stand- 
ing wave. Single frequency modulation of the intensity of 
the standing wave leads to an effective Hamiltonian for 
the center-of-mass motion given by 
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(1 -2e sin (tut + (f>)) sin 2 (kx) (1) 
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where the effective Rabi frequency is f2 e ff = SI = 

Ty/ ///sat is the resonant Rabi frequency, e is the modu- 
lation parameter, w is the modulation angular frequency, 
F is the inverse spontaneous lifetime, 5 is the detuning 
of the standing wave, t is the time, p x the momentum 
component of the atom along the standing wave, and k 
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is the wave number. Here I is the spatial-mean of the 
intensity of the unmodulated standing wave (which is 
half of the peak intensity) so O = r^/7p Ca k/2/sat where 
Isat — hcT/X 3 is the saturation intensity. A is the wave- 
length of the standing wave. 4> determines the start phase 
of the amplitude modulation. Using scaled variables Q 
the Hamiltonian is given by 



H = p 2 /2 + 2k(1 



where TL = (Ak 2 /muj 2 )H, q = 2kx 7 and p 
The driving amplitude is given by 



2esin(r-|-(/)))sm 2 (( ? /2) (2) 
(2k/muj)p x . 



2uj 2 m 



(3) 



where oj r = hk 2 /2m is the recoil frequency, r = tuj is 
the scaled time variable and Uo is the well depth. The 
commutator of scaled position and momentum is given 
by 



\p,q] = it, 



(4) 



where the scaled Planck's constant is = 8uj r /ui. For 
k = 1.2 and e = 0.20 the classical Poincare surface of sec- 
tion is shown in Fig. ^ Two symmetric regular regions 
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FIG. 1: Poincare section for a classical particle in an 
amplitude-modulated optical standing wave. Momentum and 
position (one well of the standing wave) of the particle along 
the standing wave are plotted stroboscopically with the stro- 
boscope period being equal to the modulation period. The 
central region consists of small amplitude motion. Chaos 
(dotted region) separates this region from two period- 1 re- 
gions of regular motion (represented in the Poincare section 
as sets of closed curves) located left and right of the centre 
along momentum p = . Further out in momentum are two 
stable regions of motion known as librations. At the edges 
are bands of regular motion corresponding to above barrier 
motion. It is plotted for modulation parameter e = 0.20 and 
scaled well depth « = 1.20. 

can be observed about (q — Q,p — 1) and (q — 0,p = —1). 
These regions correspond to oscillatory motion in phase 
with the amplitude modulation in each well of the stand- 
ing wave. In the experiment 0,0] atoms are loaded in a 
period- 1 region of regular motion by controlling their ini- 
tial position and momentum and by choosing the starting 



phase of the amplitude modulation appropriately. Clas- 
sically atoms should retain their momentum state when 
observed stroboscopically (time step is one modulation 
period). A distinct signature of dynamical tunnelling 
is a coherent oscillation of the stroboscopically observed 
mean momentum as shown in Fig. 0and reported in ref. 


In Sec. [H]we introduce the theoretical tools to analyze 
dynamical tunnelling by discussing Gross-Pitaevskii sim- 
ulations and the appropriate Floquet theory. We present 
a thorough analysis of the experiments from ref. 
in Sec. IIIII After showing some theoretical results for 
the experimental parameters we give a small overview 
of what to expect when some of the system parameters 
in the experiments are varied in Sec. IIVI and give some 
initial analysis. In Sec. 0we point to pathways to ana- 
lyze the quantum-classical transition for our experimen- 
tal system and give conclusions in Sec. IVII 



II. THEORETICAL ANALYSIS OF THE 
DYNAMIC EVOLUTION OF A BOSE-EINSTEIN 
CONDENSATE 

A. Dynamics using the Gross-Pitaevskii equation 

The dynamics of a Bose-Einstein condensate in a time- 
dependent potential in the mean-fieldhmit are described 
by the Gross-Pitaevskii equation 



ih 



(r, t) 
dt 



V 2 + V trap (r,i)+U(r,t) 



2m 
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+ N |*(r,i)| 2 



*(r,i) (5) 



where N is the mean number of atoms in the condensate, 
a is the scattering length with a = 2.8 nm for sodium. 
Vtrap ( r ) is the trapping potential which is turned off dur- 
ing the interaction with the standing wave and V (r, t) 
is the time dependent optical potential induced by the 
optical standing wave. The Gross-Pitaevskii equation 
is propagated in time using a standard numerical split- 
operator, fast Fourier transform method. The size of 
the spatial grid of the numerical simulation is chosen to 
contain the full spatial extent of the initial condensate 
(therefore all the populated wells of the standing wave), 
and the grid has periodic boundary conditions at each 
side (a few unpopulated wells are also included on each 
side) . 

To obtain the initial wave function a Gaussian test 
function is evolved by imaginary time evolution to con- 
verge to the ground state of the stationary Gross- 
Pitaevskii equation. Then the standing wave is turned 
on adiabatically with V (r, t) approximately having the 
form of a linear ramp. After the adiabatic turn-on, the 
condensate wave function is found to be localized at the 
bottom of each well of the standing wave. The standing 
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wave is shifted and the time-dependent potential now has 
the form 

V(r,t) = ^£(l-2esin(^ + 0))sin 2 (fca; + ^) (6) 

where ip is the phase shift which is applied to selectively 
load one region of regular motion 0] . The position repre- 
sentation of the atomic wave function (r)| 2 just before 
the modulation starts (and after the phase shift) is shown 
in Fig. El (t = OT). 

The Gross-Pitaevskii equation is used to model the 
experimental details of a Bose-Einstein condensate in an 
optical 1-D lattice. The experiment effectively consists 
of many coherent single-atom experiments. The coher- 
ence is reflected in the occurrence of diffraction peaks in 
the atomic momentum distribution (see Fig. 0J. Uti- 
lizing the Gross-Pitaevskii equation, the interaction be- 
tween these single-atom experiments is modelled by a 
classical field. Although ignoring quantum phase fluctu- 
ations in the condensate, the wave nature of atoms is still 
contained in the Gross-Pitaevskii equation and dynami- 
cal tunnelling is a quantum effect that results from the 
wave nature of the atoms. The assumption for a common 
phase for the whole condensate is well justified for the ex- 
perimental conditions as the timescales of the experiment 
and the lattice well depth are sufficiently small. It will 
be shown in the following (see Fig. |5J that the kinetic 
energy is typically of the order of 10° Hz which is much 
larger than the non-linear term in the Gross-Pitaevskii 
equation [5] which is on the order of 400 Hz. The experi- 
mental results, in particular dynamical tunnelling could 
therefore be modelled by a single particle Schrodinger 
equation in a one-dimensional single well with periodic 
boundary conditions. Nevertheless the Gross-Pitaevskii 
equation is used to model all the experimental details of 
a Bosc-Einstcin condensate in an optical lattice to guar- 
antee maximum accuracy. We will discuss and compare 
the Gross-Pitaevskii and the Floquet approaches below 
(last paragraph of section II). 

Theoretical analysis of the dynamical tunnelling exper- 
iments will be presented in this paper utilizing numeri- 
cal solutions of the Gross-Pitaevskii equation. Further- 
more we will analyze the system parameter space which 
is spanned by the scaled well depth k, the modulation 
parameter e and the scaled Planck's constant ft. In fact 
variation of the scaled Planck's constant in the simula- 
tions allows one to move the quantum system towards 
the classical limit. 



B. Floquet analysis 

The quantum dynamics of a periodically driven Hamil- 
tonian system can be described in terms of the eigen- 
states of the Floquet operator F, which evolves the sys- 
tem in time by one modulation period. In the semi- 
classical regime, the Floquet eigenstates can be associ- 
ated with regions of regular and irregular motion of the 



classical map. However when h is not sufficiently small 
compared to a typical classical action the phase-space 
representation of the Floquet eigenstates do not neces- 
sarily match with some classical (regular or irregular) 
structures^! However, initial states localized at the 
stable region around a fixed point in the Poincare section 
can be associated with superpositions of a small number 
of Floquet eigenstates. Using this state basis, one can re- 
veal the analogy of the dynamical tunnelling experiments 
and conventional tunnelling in a double well system. Two 
states of opposite parity which can be responsible for the 
observed dynamical tunnelling phenomenon are identi- 
fied. Floquet states are stroboscopic eigenstates of the 
system. Their phase space representation therefore pro- 
vides a quantum analogue to the classical stroboscopic 
phase space representation, the Poincare map. 

Only very few states are needed to describe the evo- 
lution of a wave packet that is initially strongly local- 
ized on a region of regular motion. A strongly localized 
wavepacket is used in the experiments (strongly localized 
in each well of the standing wave) making the Floquet 
basis very useful. In contrast, describing the dynamics 
in momentum or position representation requires a large 
number of states so that some of the intuitive under- 
standing which one can obtain in the Floquet basis is 
impossible to gain. For example, the tunnelling period 
can be derived from the quasi-eigenenergies of the rele- 
vant Floquet states, as will be shown below. 

For an appropriate choice of parameters the phase 
space exhibits two period-1 fixed points, which for a suit- 
able Poincare section lie on the momentum axis at ipo, 
as in reference [lj. For certain values of the scaled well 
depth k and modulation parameter e there are two dom- 
inant Floquet states \<fi±) that are localized on both fixed 
points but are distinguished by being even or odd eigen- 
states of the parity operator that changes the sign of 
momentum. A state localized on just one fixed point is 
therefore likely to have dominant support on an even or 
odd superposition of these two Floquet states: 

|^(±po))«(|0+)±|0_))/V2. (7) 

The stroboscopic evolution is described by repeated 
application of the Floquet operator. As this is a unitary 
operator, 

F\<f>±) =e ( - j2 ^±/% ± ). (8) 

<p± are the Floquet quasienergies. Thus at a time which 
is n times the period of modulation, the state initially 
localized on +po, evolves to 

\ip (n)) « U- i2 ™*+^4> + )+e(- a ™*-W\<l)-)) /y/2. 

Ignoring an overall phase and defining the separation be- 
tween Floquet quasienergies as 

A0 = 0_-0+, (10) 
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one obtains 

\1> W) « (|0+) + e (-^/*)|^_>) /V2. (n) 

At 

n = £/2A</> (12) 

periods, the state will form the anti-symmetric superpo- 
sition of Floquet states and thus is localized on the other 
fixed point at — Pq. In other words the atoms have tun- 
nelled from one of the fixed points to the other. This 
is reminiscent of barrier tunnelling between two wells, 
where a particle in one well, in a superposition of sym- 
metric and antisymmetric energy eigenstates, oscillates 
between wells with a frequency given by the energy dif- 
ference between the eigenstates. 

Tunnelling can also occur when the initial state has 
significant overlap with two non-symmetric states. For 
example, if the initial state is localized on two Floquet 
states, one localized inside the classical chaotic region 
and one inside the region of regular motion, a distinct 
oscillation in the stroboscopic evolution of the mean mo- 
mentum may be visible. The frequency of this tunnelling 
oscillation depends on the spacing of the corresponding 
quasi-eigenenergies in the Floquet spectrum. In many 
cases multiple tunnelling frequencies occur in the strobo- 
scopic evolution of the mean momentum some of which 
are due to tunnelling between non-symmetric states. 

Quantum dynamical tunnelling may be defined in that 
a particle can access a region of phase space in a way 
that is forbidden by the classical dynamics. This im- 
plies that it crosses a Kolmogorov-Arnold-Moser (KAM) 
surface [l2L IT^ | . The clearest evidence of dynamical tun- 
nelling can be obtained by choosing the scaled Planck's 
constant ft sufficiently small so that the atomic wavefunc- 
tion is much smaller than the region of regular motion 
(the size of the wave function is given by ft). Further- 
more it should be centered inside the region of regular 
motion. However, even if the wave packet is larger than 
the region of regular motion and also populates the clas- 
sical chaotic region of phase space one can still analyze 
quantum-classical correspondence and tunnelling. One 
assumes a classical probability distribution of point par- 
ticles with the same size as the quantum wave function 
and compares the classical evolution of this point parti- 
cle probability distribution with the quantum evolution 
of the wave packet. A distinct difference between the two 
evolutions results from the occurrence of tunnelling as- 
suming that the quantum evolution penetrates a KAM 
surface visibly. 



C. Comparison 

Using the Gross-Pitaevskii equation one can exactly 
simulate the experiment because the momentum or po- 
sition representation is used. Therefore the theoretical 



simulation can be directly compared with the experi- 
mental result. In contrast Floquet states do not have a 
straightforward experimental intuitive analogue. In the 
Floquet states analysis one can compare the quasienergy 
splitting between the tunnelling Floquet state with the 
experimentally measured tunnelling period. The occur- 
rence of multiple frequencies in the experimentally ob- 
served tunnelling oscillations might also be explained 
with the presence of more than two dominant Floquet 
states, the tunnelling frequencies being the energy split- 
ting between different participating Floquet states. Us- 
ing the Gross-Pitaevskii approach one can simulate the 
experiment with high precision. The same number of 
populated wells as in the experiment can be used and 
the turn-on of the standing wave can be simulated using 
an appropriate turn-on Hamiltonian. Therefore one can 
obtain the correct initial state with high precision. Fur- 
thermore the mean-field interaction is simulated which is 
not contained in the Floquet approach. 



III. NUMERICAL SIMULATIONS OF THE 
EXPERIMENTS 

A. Gross-Pitaevskii simulations 

In this section experimental results that were discussed 
in a previous paper p] are compared with numerical sim- 
ulations of the Gross-Pitaevskii (GP) equation. 

In the experiment Q] the atomic wave function was 
prepared initially to be localized around a period- 1 re- 
gion of regular motion. Figure shows the stroboscop- 
ically measured mean momentum as a function of the 
interaction time with the standing wave for modulation 
parameter e = 0.29, scaled well depth k = 1.66, mod- 
ulation frequency u>/2tt = 250 kHz and a phase shift of 
ip = 0.21 • 2tt. (a) and (b) correspond to the two differ- 
ent interaction times n + 0.25 (standing wave modula- 
tion ends at a maximum of the amplitude modulation) 
and n + 0.75 modulation periods (standing wave modu- 
lation ends at a minimum of the amplitude modulation), 
respectively. Results from the simulations (solid line, cir- 
cles) are compared with the experimental data (dashed 
line, diamonds). Dynamical tunnelling manifests itself 
as a coherent oscillation of the stroboscopically observed 
mean momentum. This occurs in contrast to the classical 
prediction in which atoms should retain their momentum 
state when observed stroboscopically (time step is one 
modulation period). There is good agreement between 
experiment and theory as far as the tunnelling period is 
concerned. However, the experimentally measured tun- 
nelling amplitude is smaller than the theoretical predic- 
tion. 

It should be noted that the theoretical simulations 
do not take account of any possible spatial and tem- 
poral variations of the scaled well depth (eg. light in- 
tensity), which could possibly lead to the observed dis- 
crepancy. It was decided to produce simulations without 
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FIG. 2: Stroboscopic mean momentum as a function of the 
interaction time with the modulated standing wave measured 
in modulation periods, n, for modulation parameter e = 0.29, 
scaled well depth k — 1.66, modulation frequency ui/2n — 250 
kHz and a phase shift ip — 0.21 ■ 2tt. (a) and (b) correspond 
to the two different interaction times n + 0.25 and n + 0.75 
modulation periods, respectively. Results from the dynamic 
evolution of the Gross-Pitaevskii equation are plotted as solid 
line (circles) and the experimental data are plotted as dashed 
line (diamonds). 



using any free parameters. However the uncertainty asso- 
ciated with the experiment would allow small variations 
of the modulation parameter e and the scaled well depth 
K. There was a 10% uncertainty in the value of the scaled 
well depth k and a 5% uncertainty in the modulation pa- 
rameter e (all reported uncertainties are 1 s.d. combined 
systematic and statistical uncertainties). Both temporal 
and spatial uncertainty during one run of the experiment 
are contained in these values as well as the systematic to- 
tal measurement uncertainty. It was verified that there 
are no important qualitative changes when varying the 
parameters in the uncertainty regime for the simulations 
presented here. However, the agreement between experi- 
ment and theory often can be optimized (not always non- 
ambiguously, meaning that it is sometimes hard to decide 
which set of parameters produces the best fit). Although 
the theoretical simulation presented does not show any 
decay in the mean momentum curve, a slight change of 
parameters inside the experimental uncertainty can lead 
to decay which is most likely caused by another dominant 
Floquet state whose presence leads to the occurrence of 
a beating of the tunnelling oscillations which appears as 
decay (and revival on longer time scales). A detailed 
analysis of the corresponding Floquet states and their 
meaning will be presented in Sec. IIII CI An example of 
how a change of parameters inside the experimental un- 
certainty in the simulations can optimize the agreement 
between theory and experiment is shown in Fig. and 
Fig. [7| Both figures are described in more detail later 
in this section. Another reason for the observed discrep- 
ancy could be the evolution of non-condensed atoms that 
is not contained in the GP approach and the interaction 
of non-condensed atoms with the condensate. With a suf- 



ficiently long adiabatic turn-on time of the far detuned 
standing wave the production of non-condensed atoms 
should be negligible. The interaction of the condensate 
with non-condensed atoms should also be negligible due 
to the low atomic density (note that in the experiments 
the condensate is expanded before the standing wave is 
turned on). However, further studies are needed to give 
an exact estimate of these effects. 

The position representation of the atomic wave func- 
tion \H! (x) | 2 is plotted stroboscopically after multiples of 
one modulation period in Fig. [2 for the same parameters 
as Fig. |21 The position of the standing wave wells is 
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FIG. 3: Position representation of the atomic wave function as 
a function of the number of modulation periods calculated us- 
ing the Gross-Pitaevskii equation. The wave packet is plotted 
stroboscopically at a phase so it is located classically approxi- 
mately at its highest point in the potential well. The position 
of the standing wave wells is also shown (dotted line). Dy- 
namical tunnelling can be observed. At t — 5T most atoms 
have tunnelled into the other period-1 region of regular mo- 
tion. The position axis is given in Fermi units (scaled by the 
mean Thomas-Fermi diameter). 

also shown (dotted line), their amplitude is given in ar- 
bitrary units. The position axis is scaled with the mean 
Thomas-Fermi diameter. The initial modulation phase is 
chosen so that the wavepacket should be located classi- 
cally approximately at its highest point of the potential 
well 0| • Choosing this stroboscopic phase the two re- 
gions of regular motion are always maximally separated 
in position space. Using this phase for the stroboscopic 
plots enables the observation of dynamical tunnelling in 
position space as the two regions of regular motion are 
located to the left and to the right of the minimum of 
the potential well, being maximally spatially separated. 
In contrast, in the experiments, tunnelling is always ob- 
served in momentum space (the standing wave is turned 
off when the regions of regular motion are at the bottom 
of the well, overlapping spatially but having oppositely 
directed momenta) as it is difficult to optically resolve 
individual wells of the standing wave. The first picture 
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in Fig. [SI (t = OT) shows the initial wave packet before 
the modulation is turned on. Subsequent pictures exem- 
plify the dynamical tunnelling process. At t = 2T, half 
the atoms have tunnelled; most of the atoms are in the 
other region of regular motion at t — 5T. The atoms 
have returned to their initial position at about t = 9T. 
The double peak structure at t = 5T and the small cen- 
tral peak at t = 2T could indicate that Floquet states 
other than the two dominant ones are also loaded which 
is likely as a relatively large Planck's constant was used 
enabling the initial wave packet to cover a substantial 
phase space area. 

Figure 01 shows simulations of the stroboscopically 
measured momentum distributions | (p) | 2 as bar graphs 
for the same parameters as Fig. |3 Corresponding exper- 
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FIG. 4: Momentum distributions as a function of the interac- 
tion time with the modulated standing wave calculated using 
numerical solutions of the Gross-Pitaevskii equation for the 
same parameters as Fig. |3] Initially atoms have mostly neg- 
ative momentum (t = 0.25T). After approximately 5 modu- 
lation periods most atoms populate a state with positive mo- 
mentum, therefore having undergone dynamical tunnelling. 
Corresponding experimental data from reference Q is also 
included as insets. 

imental data from reference [ij are also included as insets. 
The momentum distributions are plotted at 71+0.25 mod- 
ulation periods, where n is an integer. At this modulation 
phase the amplitude modulation is at its maximum and 
atoms in a period- 1 region of regular motion are classi- 
cally at the bottom of the well having maximum momen- 
tum. The two period- 1 regions of regular motion can be 
distinguished in their momentum representation at this 
phase as they have opposite momenta. At t = 0.25T the 
atoms which were located initially half way up the po- 
tential well as shown in Fig. [31 at t — 0T have "rolled" 
down the well having acquired negative momentum. Mo- 
mentum distributions for subsequent times illustrate the 
dynamical tunnelling process. At t = 5.25 modulation 
periods most atoms have reversed their momentum and 



they return to their initial momentum state at approx- 
imately t = 9.25T. The simulations arc in reasonable 
agreement with the experimentally measured data. 

Dynamical tunnelling is sensitive to the modulation 
parameter, the scaled well depth and the scaled Planck's 
constant. To illustrate this, tunnelling data along with 
the appropriate evolution of the Gross Pitaevskii equa- 
tion will be shown for another two parameter sets. Even 
though this represents only a small overview of the pa- 
rameter dependency of the tunnelling oscillations it may 
help to appreciate the variety of features in the atomic 
dynamics. Figure [S] shows the theoretical simulation and 
the experimental data for e — 0.28, k — 1.49, lu/2it = 250 
kHz, and <p = 0.22 • 27r. The mean momentum is plotted 
stroboscopically with the intensity modulation at maxi- 
mum (n + 0.25 modulation periods, n being an integer). 
The solid line (circles) is produced by a Gross-Pitaevskii 
simulation and the dashed line (diamonds) consists of ex- 
perimental data. There are approximately 3.5 tunnelling 
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FIG. 5: Mean momentum as a function of the interaction time 
with the modulated standing wave measured in modulation 
periods, n, for e = 0.28, k = 1.49, u/2tt = 250 kHz, and 
ip = 0.22 • 2-7T. The points are plotted stroboscopically with an 
interaction time of n + 0.25 modulation periods which corre- 
sponds to turning off the standing wave at maximum. Results 
from the dynamic evolution of the Gross-Pitaevskii equation 
are plotted as solid line (circles) and the experimental data is 
plotted as dashed line (diamonds). 

periods in 40 modulation periods in the theoretical curve. 
Figure shows the mean momentum as a function of 
the interaction time with the standing wave for modu- 
lation parameter e = 0.30, scaled well depth k = 1.82, 
modulation frequency lo/2tt — 222 kHz and phase shift 
ip = 0.21 ■ 27r. For these parameters the tunnelling fre- 
quency is larger than for the parameters shown in Fig. 

The simulation shows good agreement with the exper- 
iment. However, the theoretical mean momentum tun- 
nelling amplitude is larger than the one measured in the 
experiment and the theoretical tunnelling frequency for 
this set of parameters is slightly larger than the experi- 
mentally measured one. Figure \7\ illustrates that one can 
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FIG. 6: Mean momentum as a function of the interaction 
time with the modulated standing wave measured in modu- 
lation periods, n, for modulation parameter e — 0.30, scaled 
well depth k = 1.82, modulation frequency ui/2n — 222 kHz 
and phase shift <p — 0.21 • 2ir. The points are plotted stro- 
boscopically with an interaction time of n + 0.25 modulation 
periods which corresponds to turning off the standing wave at 
maximum. Results from the dynamic evolution of the Gross- 
Pitaevskii equation are plotted as solid line (circles) and the 
experimental data is plotted as dashed line (diamonds). 
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FIG. 7: Mean momentum as a function of the interaction 
time with the modulated standing wave measured in modula- 
tion periods, n, using the same parameters as Fig. |S|but the 
modulation parameter e is reduced from 0.30 to 0.28 for the 
theoretical simulation. A much better fit is obtained. 



achieve much better agreement between experiment and 
simulation if one of the parameters is varied inside the ex- 
perimental regime of uncertainty. The theoretical curve 
in this figure is obtained using the same parameters as in 
Fig. but the modulation parameter e is reduced from 
0.30 to 0.28. The tunnelling amplitude and frequency is 
now very similar to the experimental results. Note that 
the experimental data is not centered at zero momen- 
tum. This is also the case in the theoretical simulations 
for both e = 0.28 and e = 0.30. The mean momen- 
tum curve appears much more sinusoidal than the one 
for s = 0.29, k = 1.66, u/2n = 250 kHz and ip = 0.21 • 2tt 



(Fig. 0). This could imply that the initial wave function 
has support on fewer Floquet states. 



B. Stroboscopic evolution of the system energies 

Calculating the expectation values of the relevant sys- 
tem energies can give important information about the 
relevant energy scales and it might also help to obtain a 
deeper insight into the stroboscopic evolution of a Bosc- 
Einstein condensate in a periodically modulated poten- 
tial. Figure [S] shows the energy expectation values for 
the mean-field energy, the potential energy and the ki- 
netic energy given in Hz (scaled by Planck's constant). 
For comparison the stroboscopic mean momentum eX- 
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FIG. 8: Stroboscopic evolution of mean-field energy, the ki- 
netic energy and the potential energy. The energies are given 
in Hz (energy is scaled with Planck's constant). For compari- 
son the mean momentum is also shown as a function of the in- 
teraction time with the standing wave (in units of modulation 
periods). The solid (circles) and dashed (diamonds) curves 
correspond to the two different interaction times n + 0.25 and 
n + 0.75 modulation periods, respectively. The evolution is 
plotted for modulation parameter e = 0.29, scaled well depth 
k = 1.66, modulation frequency uj/2h = 250 kHz and a phase 
shift ip = 0.21 ■ 2tt. 

pectation values are also shown. The energy and mo- 
mentum expectation values are plotted stroboscopically 
and the solid (circles) and dashed (diamonds) curves cor- 
respond to the two different interaction times n + 0.25 
and n + 0.75 modulation periods, respectively. Figure |HI 
is plotted for modulation parameter e = 0.29, scaled well 
depth k = 1.66, modulation frequency u>/2tt = 250 kHz 
and a phase shift of <p = 0.21 • 2-7T which corresponds to 
Fig. [3 The mean-field energy is three orders of mag- 
nitude smaller than the potential or the kinetic energy. 
While the kinetic energy does not show a distinct oscil- 
lation, an oscillation is clearly visible in the stroboscopic 
potential energy evolution. This oscillation frequency is 
not equal to the tunnelling frequency but it is smaller as 
can be seen in Fig. [H] One obtains a period of approxi- 
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mately 8.9 modulation periods compared to a tunnelling 
period of approximately 10.0 modulation periods. Con- 
sidering the energy scales one should note that this oscil- 
lation is also clearly visible in the stroboscopic evolution 
of the total energy of the atoms. The origin of this oscil- 
lation is not known yet and will be the subject of future 
investigation. 



C. Floquet analysis for some experimental 
parameters 

Time and spatial periodicity of the Hamiltonian allow 
utilization of the Bloch and Floquet theorems 0, [w| . 
Because of the time periodicity, there still exists eigen- 
states of the evolution operator over one period (Flo- 
quet theorem). Its eigenvalues can be written in the 
form e - l27 "t>^T/k w h ere n ig ca ll ec | the quasienergy of 

the states, n is a discrete quantum number. Due to the 
spatial A/2-periodicity, in addition to n, these states are 
labeled by a continuous quantum number, the so called 
quasi-momentum $ E [— 2ir/\, 2tt/X[ (Bloch theorem). 
The quasienergy spectrum <^ n ($) is therefore made of 
bands labelled by n (see for instance Fig. I10[l . More 
precisely, the states can be written as 

\<t>n,V (r)) = e ("«*M*W*) e (-««) (r)) (13) 

where {{ipn,^ ( r ))} is now strictly periodic in space and 
time (i.e. not up to a phase). 

The evolution of the initial atomic wave function can 
be easily computed from its expansion on the \4> n ,# (0)) 
once the Floquet operator has been diagonalized. The 
iV'n.d ( T )) are the eigenstates of the modified Floquet- 
Bloch Hamiltonian 

H= (p + M) 2 /2 + 2K(l-2esinT)sin 2 ( g /2) (14) 

subjected to strictly periodic space-time boundary con- 
ditions. 

Dominant Floquet states may be determined by calcu- 
lating the inner product of the Floquet states with the 
initial atomic wavefunction. To obtain a phase space rep- 
resentation of Floquet states in momentum and position 
space one can calculate the Husimi or Q-function. It is 
defined as 

Q(q,P,r) = -L|(g + i p |$| 2 (15) 

where | q + ip) is the coherent state of a simple harmonic 
oscillator with frequency w chosen as yfH in scaled units. 
The position representation (q'\q + ip) of the coherent 
state [lj, [ljj is given by 

1/4 ( 
exp I 



(q'\q + ip) = [^ 



2A 



q 



ip| (16) 



up to an overall phase factor and where A = J . 

Floquet analysis will be shown for some of the exper- 
imental parameters which were presented in the previ- 
ous section. Figure El shows contour plots of the Husimi 




FIG. 9: Phase space representations of two Floquet states 
for 1? = that are involved in the dynamical tunnelling are 
shown in (a) and (b). The Floquet states correspond to the 
experimental parameters: modulation parameter e = 0.29, 
scaled well depth k = 1.66 and modulation frequency ui/2n = 
250 kHz, which were utilized to obtain the experimental re- 
sults shown in Fig. |5] A third state, shown in (c) also has 
significant overlap with the initial experimental state. 



functions of two Floquet states with opposite parity ((a), 
(b)) for these parameters whose presence allows dynami- 
cal tunnelling to occur. Both of them are approximately 
localized on the classical period-1 regions of regular mo- 
tion and they were selected so that the initial atomic 
wave function has significant overlap with them (26% 
and 44%, respectively). The initial experimental state 
has also significant overlap (22%) with a third state that 
is shown in Fig. Etc). The overlap is calculated using 
a coherent state that is centered on the periodic region 
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of regular motion, which is a good approximation of the 
initial experimental state. The quasi-eigenenergy spec- 
trum for this set of parameters is shown in Fig. The 




-0.5 0.5 

Bloch angle 

FIG. 10: Quasi-eigenenergy spectrum for parameters modula- 
tion parameter e = 0.29, scaled well depth n = 1.66 and mod- 
ulation frequency lu/2tt — 250 kHz. The quasi-eigenenergy is 
measured in reduced energy units and it is a function of the 
Bloch angle. The Bloch angle is equal to the quasi-momentum 
multiplied with the spatial period A/2 of the lattice. Each line 
corresponds to one Floquet state labelled by n (see Ea 113ft , 
The quasi-eigenenergy of most Floquet states is strongly de- 
pendent on the quasi-momentum but it is not the case for 
the states of the tunnelling doublet. Recall that the quasi- 
momentum i? is proportional to the phase taken by a state 
when spatially translated by A/2. In other words, fixing $ 
imposes on the states some conditions on the boundary of 
one elementary cell [q,q + A/2]. If a state is localized deep 
inside a cell, changing the boundary conditions (i.e. •&) will 
not affect the state so much and the corresponding quasiener- 
gies appear as curves that are approximately parallel to the •& 
axis. However, one should expect a strong ^-dependence for 
a state that spread over at least Aq ~ A/2. 

quasi-eigenenergies are plotted as a function of the quasi- 
momentum (Bloch angle). Each line corresponds to one 
Floquet state labelled by n (see Ea ll3l) . The arrows point 
to the two Floquet states shown in Fig. |jja) and (b) that 
correspond to the tunnelling splitting in Eq. Most 
of the other lines correspond to Floquet states which lie 
in the classical chaotic phase space region. Examples of 
phase space representations of such Floquet states can be 
found in reference Note that the spectrum is 2n1i/T 
periodic in quasienergies. For quasi-momentum $ = 
the splitting between the two states, that are shown in 
Fig. is approximately 0.08 in reduced units (energy 
in frequency units, [energy/ (27rfe)]). With a modulation 
frequency uj/2-k — 250 kHz one obtains a scaled Planck's 
constant % = 0.8, therefore the tunnelling period T tun 
which follows from Eq. El is 10 which is in good agree- 
ment with the experiment. 

The quasi-momentum plays a significant role in the 
experiments. The quasi-momentum $ is approximately 



equal to the relative velocity v between the wave packet 
(before the lattice is turned on) and the lattice, v — #/m 
[17] if the standing wave is adiabatically turned on. It was 
found in that it is of importance to populate a state 
whose quasi-momentum average is equal to zero. More- 
over, quasi-momentum spread is also of importance. It 
has been shown in that if the thermal velocity distri- 
bution is too broad, then the tunnelling oscillation disap- 
pears. As can be seen in Fig. 1101 the quasi-eigenenergies 
of the two contributing Floquet states depend on the 
quasi-momentum (or Bloch angle). The tunnelling pe- 
riod (which is inversely proportional to the separation 
between these two quasi-eigenenergies) depends on the 
quasi-momentum. Using a thermal atomic cloud one ob- 
tains a statistical ensemble of many quasi-momenta as 
they initially move in random directions with respect to 
the optical lattice. Atoms localized in individual wells 
can be described by a wave packet in the plane wave 
basis and therefore they are characterized by a super- 
position of many quasi-momenta. The resulting quasi- 
momentum width washes out the tunnelling oscillations 
(see [l8L Fig. 3]). In fact in another experiment Steck et 
al. [3J found that the amplitude of the mean momentum 
oscillations resulting from a tunnelling process between 
two librational islands of stability decreased when the ini- 
tial momentum width of the atomic cloud was increased. 



Figure shows contour plots of Husimi functions for 
the two dominant Floquet states for the modulation pa- 
rameter e = 0.30, scaled well depth k = 1.82 and mod- 
ulation frequency lo/2it = 222 kHz, which corresponds 
to experimental results shown in Fig. The states are 
selected to have maximum overlap with the initial wave 
packet (38% and 44%). In contrast to the experimen- 
tal results shown in Fig. [2] there are only two dominant 
Floquet states. The quasi-eigenenergy spectrum for this 
set of parameters (not shown) reveals a level splitting of 
approximately 0.15 in reduced units, the calculated tun- 
nelling period is 6 modulation periods which is in good 
agreement with the experiment. 



When comparing the stroboscopic evolution of the 
mean momentum shown in Fig. [5] (modulation param- 
eter e = 0.29, scaled well depth k = 1.66 and modu- 
lation frequency uj/2tt — 250 kHz) with the one shown 
in Fig. (modulation parameter e = 0.30, scaled well 
depth k = 1.82 and modulation frequency ui/2ir = 222 
kHz), one finds that it is less sinusoidal. This can be 
explained in the Floquet picture. While there are three 
dominant Floquet states for the first case (Figs. |3 
(three Floquet states with significant overlap with the 
initial experimental state), there are only two dominant 
Floquet states for the second case (Figs. El lll|l resulting 
in a more sinusoidal tunnelling oscillation. 
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FIG. 11: Phase space representations of two Floquet states 
whose presence can lead to the occurrence of dynamical tun- 
nelling. The Floquet states correspond to the experimen- 
tal parameters: modulation parameter e = 0.30, scaled well 
depth k = 1.82 and modulation frequency uj/2-k = 222 kHz, 
which were utilized to obtain the experimental results shown 
in Fig. HJ 



D. Loading analysis of the Floquet superposition 
state 

The initial atomic wave packet is localized around the 
classical period- 1 region of regular motion by inducing a 
sudden phase shift to the standing wave. This enables 
the observation of dynamical tunnelling. In the Floquet 
picture the observation of dynamical tunnelling requires 
that the initial state has support on only a few domi- 
nant Floquet states, preferably populating only two with 
a phase space structure as shown in Fig. [5J Optimizing 
the overlap of the initial state with these Floquet "tun- 
nelling" states should maximize the observed tunnelling 
amplitude. Here we carry out analysis confirming this 
prediction. This corresponds to optimizing the overlap of 
the initial experimental state with the period- 1 regions 
of regular motion. Figure IT2*1 shows the mean atomic mo- 
mentum as a function of the interaction time with the 
standing wave which is plotted for a range of the initial 
phase shift ip of the standing wave. The solid line (cir- 
cles) corresponds to an interaction time with the stand- 
ing wave of n + 0.25 modulation periods and the dashed 
line (diamonds) corresponds to an interaction time of 
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FIG. 12: Mean atomic momentum as a function of the in- 
teraction time with the standing wave plotted for a range of 
the initial phase shift ip of the standing wave. The solid line 
(circles) corresponds to an interaction time with the stand- 
ing wave of n + 0.25 modulation periods and the dashed line 
(diamonds) corresponds to an interaction time of n + 0.75 
modulation periods. 



n + 0.75 modulation periods. The simulations are made 
using the Gross-Pitaevskii equation for modulation pa- 
rameter e = 0.29, scaled well depth k = 1.66, modu- 
lation frequency lo/2tt — 250 kHz which corresponds to 
Fig. [21 A phase shift ip = corresponds to localizing 
the wave packet exactly at the bottom of the well and 
(p = tt corresponds to localizing it exactly at the maxi- 
mum of the standing wave well. Symmetry dictates that 
no tunnelling oscillation can occur for these two loading 
phases. This is also shown in Fig. ^] (the mean mo- 
mentum curve for tp = tt is not shown but it is the same 
as for ip — 0). The amplitude of the tunnelling oscil- 
lations changes strongly when the initial phase shift ip 
of the standing wave is changed. The best overlap with 
the tunnelling Floquet states is obtained for the phase 
shift ip somewhere in between 0.25 x 2-7T and 0.30 x 2-7T 
which corresponds to placing the wave packet half way 
up the standing wave well. This result is in good agree- 
ment with the structure of the tunnelling Floquet states 
as shown in Fig. |5J When changing ip there is no change 
in the observed tunnelling period. This is to be expected 
as mainly the loading efficiency for the "tunnelling" Flo- 
quet states varies when ip is varied. It should be noted 
that the simulations are carried out for a relatively large 
scaled Planck's constant (fe = 0.8) which means that the 
wave packet size is rather large compared to character- 
istic classical phase space features like the period- 1 re- 
gions of regular motion. This analysis shows that the 
tunnelling amplitude sensitively depends on the loading 
efficiency of the tunnelling Floquet states and that there 
is a smooth dependency of the tunnelling amplitude on 
this loading parameter. 
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IV. PARAMETER DEPENDENCY OF THE 
TUNNELLING OSCILLATIONS 

The scaled parameter space for the dynamics of the 
system is given by the scaled well depth k and the mod- 
ulation parameter e. Both parameters will significantly 
change the structure and number of the contributing Flo- 
quet states. It has been found that a strong sensitiv- 
ity of the tunnelling frequency on the system parame- 
ters is a signature of chaos-assisted tunnelling [Til Hsj 
where a third state associated with the classical chaotic 
region interacts with the tunnelling Floquet states. A 
Floquet state that is localized inside a region of regular 
motion that surrounds another resonance can also inter- 
act with the tunnelling doublet (two tunnelling Floquet 
states), this phenomenon is known as resonance-assisted 
tunnelling 19] . 

Comprehensively exploring the parameter space and 
its associated phenomena is out of the scope of this arti- 
cle. Instead an analysis associated with our experiments 
will be presented here showing one scan of the scaled well 
depth k and another scan of of the modulation param- 
eter e around the experimental parameter regime. The 
results are shown in the form of plots of the mean mo- 
mentum as a function of the interaction time with the 
modulated standing wave. The solid line (circles) cor- 
responds to an interaction time with the standing wave 
of n + 0.25 modulation periods and the dashed line (di- 
amonds) corresponds to an interaction time of n + 0.75 
modulation periods. 

Figure IT51 shows the scaled well depth n being varied 
from 1.10 to 1.75. The other parameters are held con- 
stant (modulation parameter e = 0.29, modulation fre- 
quency oj/2tt — 250 kHz and phase shift ip — 0.21 • 2tt). 
The momentum distribution evolution is shown to illus- 
trate intricate changes in the tunnelling dynamics when 
the parameters are varied. Both the tunnelling frequency 
spectrum and the tunnelling amplitude are strongly de- 
pendent on k. Often one can see more than just one dom- 
inant tunnelling frequency. For example, for n = 1.35 the 
two dominant tunnelling frequencies which contribute to 
the tunnelling oscillation have a period of approximately 
3.9 modulation periods and approximately 34 modula- 
tion periods. In the interval of approximately n = 1.5 
and k = 1.8 the tunnelling oscillations have a more si- 
nusoidal shape indicating the presence of only approxi- 
mately two dominant tunnelling states. In this interval 
the tunnelling frequency is peaked at k = 1.55 with a tun- 
nelling period of approximately 8.9 modulation periods. 
In order to analyze further the behavior of the tunnelling 
frequency, we show the low frequency (diamonds) and the 
high frequency component (squares) of the tunnelling fre- 
quency as a function of the scaled well depth k in Fig. 1141 
The high frequency component is plotted only for values 
of k where it is visible in the stroboscopic momentum 
evolution shown in Fig. 1131 The occurrence of multiple 
tunnelling frequencies results from the presence of more 
than two dominant Floquet states. The minimum of the 
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FIG. 13: Mean momentum as a function of the interaction 
time with the modulated standing wave for different values 
of the scaled well depth k (modulation parameter e = 0.29, 
modulation frequency lu/2tt = 250 kHz and phase shift ip = 
0.21 -2tt). The solid line (circles) corresponds to an interaction 
time with the standing wave of n + 0.25 modulation periods 
and the dashed line (diamonds) corresponds to an interaction 
time of n + 0.75 modulation periods. 
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FIG. 14: Tunnelling frequency as a function of the scaled 
well depth k (modulation parameter e = 0.29, modulation 
frequency ui/2n = 250 kHz and phase shift ip — 0.21 ■ 2n) 
taken from Fig. 1131 Diamonds and squares denote low and 
high frequency components, respectively. 



12 




FIG. 15: Floquet spectrum as a function of the scaled well depth k (modulation parameter e = 0.29, modulation frequency 
uj/2-k — 250 kHz and phase shift ip — 0.21 • 2tt. The Floquet states with maximum overlap are marked black bullets, white 
bullets are used for second most overlap and white squared bullets show third most overlap. Phase space representations 
(Husimi functions) of some of the these Floquet states for different values of k are also shown. 



low frequency component of the tunnelling frequency is 
due the the level crossing of two contributing Floquet 
states (see squares and filled circles at k = 1.25 in Fig. 
I15JI . The error bars result from the readout uncertainty of 
the tunnelling period from the simulations. This analy- 
sis reveals some of the intricate features of the tunnelling 
dynamics. Instead of a smooth parameter dependency a 
distinct rise and fall (in the low frequency component) 
in the tunnelling frequency versus scaled well depth ap- 
pears. The tunnelling frequency minimum is centered 
at K w 1.3. Note that a rise and fall in the tunnelling 
frequency is often understood as a signature of chaos- 
assisted tunnelling. However, it is not a sufficient crite- 
rion for chaos-assisted tunnelling. One needs to choose 
an approximately ten times smaller scaled Planck's con- 
stant to use the terminology of chaos-assisted tunnelling 
[l8[ . The size of the Floquet states is given by the scaled 
Planck's constant. If the states are much larger than 
phase space features like regions of regular motion, then 
it is impossible to make a classification of Floquet states 
as chaotic or regular, which is needed for chaos-assisted 
tunnelling. 

Another interesting feature of dynamical tunnelling 
can be derived from Fig. It shows that for a cer- 
tain parameter regime the tunnelling frequency decreases 
with decreasing scaled well depth. This is the contrary 



of what one would expect for spatial energy barrier tun- 
nelling. This feature has been also observed in recent ex- 
periments |4| . Floquet spectra provide alternative means 
of analyzing the tunnelling dynamics. The dependence of 
the tunnelling frequency on the scaled well depth k can 
be understood using the appropriate Floquet spectrum. 
Figure El shows the quasi-eigenenergies of different Flo- 
quet states as a function of the scaled well depth k. The 
Floquet states with maximum overlap are marked with 
bullets. Figure ^3 also shows phase space representa- 
tions (Husimi functions) of some of these Floquet states 
for different values of k. The shape and structure of 
these Floquet states depend on the value of the scaled 
well depth. In fact Floquet states can undergo bifurca- 
tions. This may be seen as the quantum analogue of 
classical phase bifurcation. Classical phase space bifur- 
cations have been reported in ref. |2fJ. In this case the 
shapes of the Floquet states change in such a way that 
different Floquet states have non-negligible overlap with 
the initial experimental state as the scaled well depth k is 
varied. The separation between the quasi-eigenenergies 
of the two states with maximum overlap will determine 
a dominant tunnelling frequency. Note that often more 
than two states have relevant overlap with the initial ex- 
perimental states which leads to the occurrence of mul- 
tiple tunnelling frequencies. 
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Figure ^] shows effects of a variation of the modula- 
tion parameter e. To obtain these simulations all other 




should enhance our understanding of nonlinear dynami- 
cal systems and provide insight into their quantum and 
classical origin. It is out of the scope of this paper to 
present anything more than a short analysis relevant to 
the experimental results. The quantum-classical border- 
land is analyzed by considering the mean momentum as 
a function of the interaction time with the modulated 
standing wave for different values of the scaled Planck's 
constant ft. Figure IT7I shows results for modulation pa- 
rameter e — 0.29, scaled well depth n = 1.66 and phase 
shift ip — 0.21 • 2tt which corresponds to experimental re- 
sults shown in Fig. [21 The scaled Planck's constant ft is 
varied by adjusting the modulation frequency lojl'K and 
leaving the scaled well depth k, the modulation parame- 
ter e and the initial phase shift ip constant. Results are 
shown for the scaled Planck's constant ft ranging from 
0.40 (w/2tt = 500 kHz) to 1.33 (w/2tt = 150 kHz). The 
error bars result from the readout uncertainty of the tun- 
nelling period from the simulations. As shown in previous 



FIG. 16: Tunnelling frequency as a function of the modu- 
lation parameter e (scaled well depth k = 1.66, modulation 
frequency ui/2ti — 250 kHz and phase shift <p — 0.21 • 2tt). 
The error bars result from the readout uncertainty of the tun- 
nelling oscillations from the simulations. The high frequency 
component that is present in the data for 0.00 < e < 0.18 is 
not shown. 

parameters are held constant (scaled well depth k = 1.66, 
modulation frequency lu/2it = 250 kHz and phase shift 
ip = 0.21 ■ 27r). For smaller values of e the correspond- 
ing classical phase space is mainly regular. There is no 
distinct oscillation at e = 0.10 which is centered at zero 
momentum. Distinct tunnelling oscillations start to oc- 
cur at e — 0.14 with a gradually increasing tunnelling 
frequency. While there is a tunnelling period of approxi- 
mately 67 modulation periods at e = 0.14, the tunnelling 
period is only approximately 7.7 modulation periods at 
e = 0.38. For larger values of the modulation parameter 
e the oscillations become less sinusoidal indicating the 
presence of an increasing number of dominant Floquet 
states. The error bars result from the readout uncer- 
tainty of the tunnelling period from the simulations. 
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FIG. 17: Dominant tunnelling frequency as a function of 
the modulation frequency uj/2h for modulation parameter 
e = 0.29, scaled well depth k = 1.66 and phase shift ip — 
0.21 • 2n which corresponds to experimental results shown 
in Fig. [5] As expected, the overall tendency consists of a 
decreasing of the tunnelling frequency as the modulation fre- 
quency increases and therefore as the scaled Planck's con- 
stant decreases, however the transition is not as smooth as 
one might expect. 



V. MOVING THE QUANTUM SYSTEM 
TOWARDS THE CLASSICAL LIMIT 

A fundamental strength of the experiments which are 
discussed here is that they are capable of exploring the 
transition of the quantum system towards the classical 
limit by decreasing the scaled Planck's constant ft. 

Here a quantum system with mixed phase space ex- 
hibiting classically chaotic and regular regions of mo- 
tion is moved towards the classical limit. By adjust- 
ing the scaled Planck's constant ft of the system, the 
wave and particle character of the atoms can be probed 
although some experimental and numerical restrictions 
limit the extent of this quantum-classical probe. This 



work [2l|. the momentum of the regions of regular mo- 
tion is proportional to the modulation frequency. Note 
that this is a purely classical feature that would disap- 
pear if the scaled momentum would be plotted instead 
of the real momentum. Intuitively one would expect the 
tunnelling frequency to decrease as the system becomes 
more classical. However, a different quite surprising phe- 
nomenon occurs. Considering the modulation frequency 
interval between 225 kHz and 325 kHz the tunnelling fre- 
quency gradually decreases as the system becomes more 
classical. The same applies for the interval between 150 
kHz to 175 kHz and the interval from 425 kHz to 475 kHz. 
Qualitative changes occur at approximately 200 kHz, 350 
kHz (375 kHz) and 500 kHz where the tunnelling oscilla- 
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tion shows first multi- frequency contributions and is then 
significantly increased as the effective Planck's constant 
is further decreased. A more detailed analysis will be 
subject of future study. The three intervals where the 
tunnelling frequency decreases appear as three arms in 
the graph. At w/2tt = 200 kHz, 350 kHz and 375 kHz 
the system cannot be well described using one dominant 
tunnelling frequency, therefore no data is shown at these 
frequencies. 

The results of this first analysis indicate that the tran- 
sition from quantum to classical physics contains many 
fascinating details to be explored. The results for the 
driven pendulum in atom optics shows that the transi- 
tion from quantum to classical dynamics is not smooth. 



K. 

However this cannot be interpreted as a signature of 
chaos-assisted tunnelling essentially because the states 
cannot be clearly identified with chaotic or regular re- 
gions of the classical phase space, therefore the notion 
of chaos-assisted tunnelling is difficult to use in this con- 
text. We note that it is important to decrease the scaled 
Planck's constant by at least an order of magnitude 0] 
for an observation of chaos-assisted tunnelling to be non- 
ambiguous. We have simulated when the system is moved 
towards the classical limit. A bifurcation-like behavior 
results which can be understood in terms of quantum 
bifurcations of the contributing Floquet states. 



VI. CONCLUSION 

Recent dynamical tunnelling experiments were ana- 
lyzed using simulations of the Gross-Pitaevskii equation 
and the corresponding Floquet theory. The main fea- 
tures of the experiments can be explained by a two or 
three state framework that is provided by Floquet the- 
ory. We have identified the relevant Floquet states and 
shown their Husimi functions. Note that the mean-field 
interaction was negligible in the experiments. Tunnelling 
period and amplitude are in good agreement with GP 
simulation and Floquet theory. 

An analysis of the parameter space has shown that 
tunnelling frequency is strongly dependent on the system 
parameters scaled well depth k, modulation parameter e 
and the scaled Planck's constant 1t. While we find an 
approximately linear dependence of the tunnelling fre- 
quency on the modulation parameter e for the set of ex- 
perimental parameters, there is a distinct spike in the 
tunnelling frequency as a function of the scaled well depth 
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